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LIDAR  Remote  sensing  of  Sound  Velocity  in  the  Ocean 

_ Technical  report  (RF-3249) 


1.  Introduction: 


Both  Raman  and  Brillouin  scattering  can  be  considered  for  sound  velocity  measurements  in 
the  ocean  using  LIDAR.  Previous  experimental  efforts  have  concentrated  on  Raman  scattering.  The 
broad  linewidth  of  Raman  radiation  (about  15  nm)  does  not  permit  good  background  noise  and 
fluorescence  rejection,  and  measurements  have  not  shown  the  required  accuracy.  While  they  have 
similar  cross  sections,  1.5  10*4  m-l  sr^  for  Raman,  2.4  1(H  m-i  sr^  for  Brillouin,  the  Brillouin 
lines  have  the  advantage  that  their  linewidth  is  only  0.5  GHz.  A  proper  design  of  the  receiver  makes 
a  system  based  on  this  scattering  process  insensitive  to  background  radiation  and  provide  accurate 
measurements.  Typical  Brillouin  shifts  are  between  7  and  8  GHz,  and  a  sensitivity  of  the  order  of  a 
MHz  is  needed  (i.e.  1  part  in  500  of  Brillouin  line). 

To  demonstrate  a  Brillouin  LIDAR  system,  the  first  year's  work  has  been  in  two  directions: 
(1)  analysis  of  the  limitation  on  the  accuracy  of  the  sound  velocity  retrieval  using  Brillouin 
scattering,  and  (2)  development  of  the  Brillouin  LIDAR  receiver  using  the  edges  of  Br2  and  I2 
absorption  lines.  The  progresses  for  each  part  are  described  in  Sect  2  and  3  respectively. 

2.  Sound  velocity  determination 


As  a  baseline  for  the  discussion,  we  consider  an  airborne  LIDAR,  operating  at  532  nm,  with  a 
energy  of  5{X)  mJ/pulse,  single  pulse  measrjement,  and  a  depth  resolution  of  1.1  m.  Accuracy  of  the 
measurement  could  be  increased  by  averaging  over  several  shots,  decreasing  the  depth  resolution, 
or  decreasing  the  LIDAR  total  range  by  using  the  system  on  a  boat,  submarine,  or  buoy.  Depending 
on  the  depth  and  the  hydrosol  content  of  the  ocean  water,  the  standard  deviation  Avb  on  the 
Brillouin  shift  measurement  range  between  0.5  and  4  MHz  [Hickman  et.  al.,  1991]. 

We  have  made  the  first  rigorous  analysis  of  the  accuracy  of  sound  velocity  measurements  using  a 
Brillouin  LIDAR.  Previous  studies  were  very  incomplete  and  frequently  misleading  (Guagliardo  et. 
al.,  1980;  Hirschberg  and  Dufilho,1984;  Leonard  and  Sweeney  1988;  and  Hickman  et.  al.,  1991).  To 
determine  the  sound  velocity  from  the  Brillouin  shift  measurement,  there  are  five  parameters 
(Salinity,  Temperature,  Sound  velocity.  Refractive  index,  and  Brillouin  shift);  past  studies  identify, 
in  the  best  case,  only  four  relationships  among  them.  To  solve  the  system  of  equations,  by  providing 
a  fifth  relationship,  we  have  analyzed  the  historical  salinity  measurements  provided  by  the  data  set  of 
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the  Oceanographic  Station  time  Series  compiled  by  the  national  Oceanographic  Data  Center 
(NODC).  More  than  75,000  data  have  been  analyzed  to  determine  the  knowledge  in  the  salinity 
mean  value  and  standard  deviation  for  a  given  location  and  for  a  given  depth.  The  important  result 
about  that  is  that  for  a  given  depth  and  location,  the  salinity  standard  deviation  is  small  and  will 
enable  us  to  retrieve  sound  velocity  with  an  accuracy  better  than  1  m/s  for  all  locations.  Such  an 
accuracy  is  sufficient  for  numerous  applications.  A  summary  of  our  results  for  sound  velocity  error 
as  a  function  of  the  salinity  standard  deviation  AS  for  various  measurements  of  the  Brillouin  shift 
accuracy  Avb  is  shown  in  Fig.  1;  a  salinity  standard  deviation  histogram  is  superimposed  on  it.  A 
paper  describing  these  results  has  been  submitted  to  Applied  Optics.  A  copy  of  it  is  attached  to  this 
report. 
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Salinity  standard  deviation  AS 

Fig  1.  Sound  velocity  standard  deviation  as  a  function  of  the  Salinity  standard  deviation 
superposed  with  the  salinity  standard  deviation  histogram,  computed  from  Ref.  7. 


The  standard  deviation  of  the  Salinity-Temperature  (S-T)  relationship  was  also  examined  using 
the  data  set  of  the  NODC .  Some  typical  plots  of  the  salinity  (in  %o)  as  a  function  of  the  temperature 
(in  *C)  are  shown  in  Fig.  2  and  3.  The  mean  value  of  the  salinity  is  represented  by  a  solid  line,  and 
the  standard  deviation  with  dotted  lines.  The  diamond-shaped  points  represent  the  measured  pairs 
(T,S).  As  can  be  seen  in  Fig.  2,  we  have  found  a  clear  correlation  between  salinity  and  temperature 
for  a  few  locations,  generally  for  large  depths.  Nevertheless,  the  general  case,  in  particular  for  the 
upper-ocean  mixed  layers,  is  almost  always  represented  by  Fig.  3  which  shows  no  correlation.  For 
our  application,  this  means  that  for  some  locations  only,  the  precision  of  the  measurement  could  be 
increased  by  considering  the  S-T  relationship. 
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Fig  2  Mean  value  (solid  line)  and  standard  deviation  (doted  lines)  of  the  salinity,  and  pairs  (T,S) 
confuted  from  Ref.  7  (diamond-shaped  points)  for  the  location  37.85N,  70.82E.  at  a  depth  of 
100  mNotice  the  correlation  between  the  salinity  and  the  temperature. 

RWO.Oe:  002$  [m]  42.88N  ie.14W 


Fig  3  Mean  value  (solid  line)  and  standard  deviation  (doted  lines)  of  the  salinity,  and  pairs  (T JS) 
computed  from  Ref.  7  (diamond-shaped  points)  for  the  location  62.00N,  3327E.  at  a  depth  of 
25  m.  There  is  no  correlation  between  the  salinity  and  the  temperature. 
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3.  Brillouin  shift  measurement 

The  typical  Brillouin  shift  sensitivity  to  be  expected  is  of  the  order  of  a  MHz.  For  LBDAR 
measurements  the  anticipated  parameters  are  a  depth  resolution  of  1.1  m  (which  requires  a  10  ns 
laser  pulse),  and  a  0.1  GHz  data  acquisition  rate.  The  technical  requirements  are:  (1)  a  stabilized 
Fourier  transform  limited  bandwidth  pulse  laser;  and  (2)  a  non  scanning  spectroscopic  method  for 
Brillouin  shift  discrimination. 


In  preliminary  work  to  be  used  as  a  reference  for  our  work  we  have  used  a  Fabry-Perot 
interferometer  to  make  the  first  experimental  frequency  resolved  measurements  of  the  Brillouin 
lines  with  a  pulsed  laser  at  high  resolution.  Fig.  4  shows  a  measured  frequency  spectrum, 
superposed  on  the  theoretical  one  for  a  water  temperature  of  29*C  (corresponding  to  a  Brillouin 
shift  of  7.55  GHz).  Only  one  free  spectral  range  of  the  Fabry-Perot  is  presented  in  the  figure.  The 
experimental  points  are  raw  data  from  a  single  scan;  no  further  processing  has  been  performed. 
There  is  very  good  agreement  with  theory.  The  strong  central  peak  is  due  to  particulate 
contamination  of  the  water.  Prior  to  these  results,  no  frequency  resolved  measurements  using  a 
pulsed  laser  had  been  made.  The  only  available  previous  results  did  not  demonstrate  sufficient 
frequency  resolution  due  to  laser  linewidth  and  probably  poor  collimation  through  the  Fabry-Perot 
frequency  analyzer.  [Guagliardo  et.  al.,  1980]. 
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These  measurements  and  the  experimental  set  up  have  been  presented  in  the  Ocean  Optics 
Xin  conference  in  Halifax,  and  have  been  published  in  the  SPEE  conference  proceedings.  A  copy  is 
attached  to  this  report.  These  results  using  a  scanning  measurement  technique  demonstrate  that  the 
stability  of  the  laser  is  not  a  limiting  factor  for  the  measurement  of  the  Brillouin  frequency  shift 
using  a  LIDAR  instrument,  the  results  also  provide  an  essential  tool  for  the  next  tasks  of  the  project. 
For  application  to  Brillouin  LDDAR  measurements,  we  require  a  better  accuracy  (the  scanning 
Fabry-Perot  technique  will  not  provide  a  measurement  accuracy  of  1/5(X)  the  Brillouin  linewidth) 
and  a  non  scanning  method  (due  to  LIDAR  timing  requirement). 

The  edge  technique  used  in  conjunction  with  Br2  and  I2  absorption  cells  will  serve  these 
requirements.  Explicitly,  the  frequency  of  the  laser  transmitter  is  chosen  so  that  the  edges  of  two 
different  molecular  absorption  lines  (with  opposite  slopes)  are  located  at  the  centers  of  the  Brillouin 
lines.  Because  of  the  steep  slope  of  the  absorption  edges,  small  changes  in  the  Brillouin  shift 
produce  large  changes  in  transmission.  A  more  precise  description  of  the  experimental  set  up  is 
given  in  the  proceeding  of  the  Halifax  conference  attached  to  this  report.  For  this  spectroscopic 
purpose,  we  succeeded  in  obtaining  second  harmonic  generation  (1064  nm  to  532  nm)  for  measuring 
the  I2  and  B12  spectra  and  for  locking  the  frequency  of  the  seed  laser  used  to  injection  lock  the 
pulsed  Nd:Yag  laser.  A  nice  feature  is  that  our  new  method  does  not  require  stabilization  of  the 
doubling  crystal  temperature,  contrary  to  past  approaches  which  required  stabilization  of  the  crystal 
to  a  few  millidegrees.  The  spectroscopic  set  up  is  now  ready.  The  Br2  cells  should  be  ready  for  the 
end  of  the  year,  and  the  *29i2  ceW  should  arrive  soon  from  China. 

3.  Summary 

The  project  is  progressing  very  well,  and  the  first  sound-velocity  measurement  using  our  non¬ 
scanning  technique  should  be  performed  by  the  summer  1997.  No  major  problems  have  been 
encountered  during  the  first  year  of  this  project.  Several  improvements  of  the  experimental  set  up, 
and  new  possibilities  are  being  considered: 

-  Faraday  Anomalous  Dispersion  Optical  Filters  (FADOF)  using  an  excited  state  of  Potassium 
could  be  also  used  as  edges  filters  and  their  use  is  under  study.  It  will  have  the  advantage  of 
rejecting  background  noise  and  having  a  very  high  transmission. 

-  The  width  of  the  Brillouin  lines  is  due  to  acoustic  damping  (absorption).  By  measuring  this 
linewidth  we  can  obtain  data  on  the  loss  mechanisms  -  specifically  bulk  viscosity  and  shear 
viscosity.  There  is  presently  no  experimental  data  for  water  on  the  former.  An  experimental  set 
up  similar  to  the  one  for  Brillouin  frequency  shift  measurement  could  serve  this  measurement. 
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Since  this  project  is  so  closely  involved  with  physical  oceanography,  the  graduate  student  in  this 

project,  Jeffery  Katz,  has  enrolled  in  a  course  in  graduate  level  physical  oceanography  this 

semester. 
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Laboratory  development  of  a  LIDAR  for  measurement  of  sound  velocity  in  the 

ocean  using  Brillouin  scattering 
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ABSTRACT 

There  have  been  several  studies  of  the  potential  accuracy  of  LIDAR  measurements  of  sound  velocity  in  the  ocean  by 
measuring  the  spectral  shift  of  the  backscattered  Brillouin  lines.  However,  due  to  technical  limitations,  such  systems  have 
not  previously  been  experimentally  demonstrated.  Measurement  of  the  Brillouin  shift  as  a  function  of  depth  in  the  ocean 
requires  a  stabilized,  narrow  linewidth,  pulsed  laser,  and  a  high-resolution  spectroscopic  technique.  We  have  used  a 
scanning  Fabry-Perot  to  obtain  the  first  frequency  resolved  measurements  of  Brillouin  scattering  in  water  using  a  pulsed 
laser,  these  results  will  be  presented  h^.  But  for  the  practical  application  of  measuring  Brillouin  shifts  as  a  function  of 
depth  in  the  ocean  a  non-scanning  spectroscopic  technique  is  required  to  measure  the  small  frequency  shifts;  the  edge 
technique  meets  this  requirement.  Using  it  in  conjunction  with  the  edges  of  absorption  lines  in  the  molecular  spectra  of  I2 
and/or  Br2.  avoids  the  limitations  associated  with  use  of  a  Fabry-Perot  etalon;  specifically,  its  small  solid  angle  of 
acceptance  and  its  vulnerability  in  noisy  environments.  This  new  approach  will  be  briefly  described. 

Keywords:  Brillouin  LIDAR  remote  sensing,  sound  velocity,  ocean,  edge  technique,  molecular  absorption  lines. 

1.  INTRODUCTION 

Several  studies  have  discussed  the  importance  of  remote  measurements  of  sound  velocity  in  the  upper-ocean  mixed 
layer. Real-time  range-resolved  monitoring  provides  a  crucial  input  to  the  understanding  of  the  physictd  behavior  of  the 
ocean.  A  Brillouin  LIDAR  is  a  promising  tool  for  such  remote  measurements  of  the  sound  velocity  as  a  function  of  depth. 
Brillouin  scattering  produces  two  inelastically  scattered  Lorentzians  centered  symmetrically  with  respect  to  the  laser  line; 
they  have  a  FWHM  of  about  0.5  GHz.  The  so-called  Brillouin  shift  vb,  that  is  to  say  the  frequency  shift  between  the  laser 
line  and  each  of  the  Brillouin  lines,  is  typically  between  7  and  8  GHz,  and  is  proportional  to  the  sound  velocity.  In  pure 
water,  the  scattoing  spectrum  is  almost  solely  composed  this  doublet  An  elastically  scattered  central  line  (also  called  the 
unshifted  line,  or  improperly  the  Rayleigh  line)^  appears  when  the  water  is  contaminated  by  hydrosols. 

The  limits  in  the  accuracy  of  sound  velocity  and  temperature  profiles  in  the  ocean  from  LIDAR  Brillouin  shift 
measurements  have  been  discussed  in  detail  by  Fry,  et  al?  Briefly,  for  a  given  incident  laser  wavelength  X  the  Brillouin 
shift  Vb  is  proportional  to  the  refractive  index  of  water  n  and  the  sound  velocity  v^ .  The  refractive  index  itself  is  function 
of  the  salinity  S  and  temp^ture  T.  In  addition  to  these  two  variables,  the  sound  velocity  also  has  a  weak  dependence  on 
pressure,  but  for  present  purposes,  this  dependence  will  be  neglected;  the  generalization  to  include  pressure  dependence  is 
straightforward.  The  Brillouin  shift  is  then  expressed  as: 

v,(S,T)  =  T)sin(9/2),  (1) 

where  6  is  the  scattmng  angle  (sin(0/2>sl  for  180*  backscattering).  Clearly,  th^e  are  five  parameters  (vg,  v,,  n,  S,  and  T), 
but  ihere  are  only  three  relations  between  them  and  only  one  measurement  (vg).  The  three  relations  are  Eq.  1  and  empirical 
equations  for  n(S,  T)  and  vs(S,  T).  Thus  one  additional  piece  of  information  is  required.  For  jvesent  purposes,  it  is 
provided  by  the  means  and  standard  deviations  of  the  extensive  set  of  salinity  measurem^ts  (from  1900  to  1990)  compiled 
by  the  National  Oceanic  and  Atmospheric  Administration.®  With  these  salinity  estimates  based  on  historical  data,  i.e. 
without  actual  measurements  of  salinity,  the  limits  on  the  sound  velocity  accuracy  can  be  better  than  0.3  m/s. 

Our  primary  interest  is  in  developing  a  LIDAR  for  range-resolved  measurements  of  sound  velocity  in  the  ocean  via  the 
measurem^t  of  the  Brillouin  shift.  Although  the  idea  of  a  Brillouin  LIDAR  is  not  new,^'^  such  systems  have  not  yet  been 
experimentally  demonstrated  due  to  technical  limitations.  Measurements  of  the  Brillouin  shift  as  a  function  of  depth  in  the 
ocean  requires  both  a  stabilized  narrow  linewidth  pulsed  laser,  and  a  high-resolution  non-scanning  spectroscqpic  technique. 


Sound  velocity  measurements  with  depth  resolution  Sz  of  1  meter  require  a  Brillouin  frequency  shift  measurement  for  each 
successive  10  ns  (=2n8zfc)  interval  of  the  LIDAR  return.  In  addition,  a  precision  of  1  m/s  for  the  sound  velocity  requires 
the  measurement  of  Brillouin  shifts  with  an  accuracy  better  than  about  4  MHz.^  This  small  value  compared  to  the 
linewidth  of  the  Brillouin  line  (500  MHz)  does  not  allow  the  use  of  a  conventional  spectrometer.  In  this  p^r  we  discuss 
our  recent  progress  in  development  of  a  Brillouin  LIDAR. 

The  first  laboratory  measurements  of  frequency  resolved  Brillouin  lines  using  a  pulsed  laser  and  a  scanning  Fabry-Perot 
technique  for  frequency  discrimination  are  presented  in  Sect.  2.  Previous  pulsed  laser  measurements^  were  not  completely 
resolved,  but  our  results  clearly  demonstrate  this  resolution  with  pulsed  lasers.  However,  the  scanning  Fabry-Perot 
measurement  technique  used  to  obtain  this  data  is  not  accurate  and  fast  enough  to  be  used  directly  for  LIDAR 
measurements  of  the  Brillouin  shift.  Nevertheless,  these  results  and  the  measurement  apparatus  will  provide  a  valuable 
reference  for  the  initial  data  obtained  with  our  new  non-scanning  spectroscopic  technique. 

This  new  technique  is  a  version  of  the  edge  technique,  and  is  a  powerful  non-scanning  spectroscopic  technique  for  the 
measurement  of  small  frequency  shifts.^  It  has  numerous  applications,  and  has  been  successfully  used  in  a  DopplCT  LIDAR 
for  atmospheric  wind  velocity  measurements.^®  From  the  standpoint  of  an  oceanic  Brillouin  LIDAR,  the  main  limitations 
of  previous  implementations  of  this  method  are  the  small  acceptance  angle  of  the  Fabry-Perot  etalons  used  for  the 
frequency  discrimination,  and  the  difficulties  associated  with  operating  these  interferometers  in  noisy  vibrational  and 
acoustic  environments.  We  are  implementing  a  version  of  the  edge  technique  in  which  the  frequency  discriminants  are  the 
edges  of  absorption  lines  in  the  molecular  spectra  of  I2  and/or  Br2.  Explicitly,  the  frequency  of  ^e  laser  transmitter  is 
chosen  so  that  the  edges  of  two  different  molecular  absorption  lines  (with  opposite  slopes)  are  located  at  the  centos  of  the 
Brillouin  lines.  Because  of  the  steep  slope  of  the  absorption  edges,  small  changes  in  the  Brillouin  shift  produce  large 
changes  in  transmission.  The  details  of  the  experimental  set  up  are  described  in  Sect  3. 

2.  FREQUENCY  RESOLVED  MEASUREMENTS  OF  BRILLOUIN  LINES 

The  theoretical  model  for  Brillouin  scattering  in  a  liquid  was  developed  by  Mountain^  wd  Fabelinskii.^^  Early  attempts 
to  examine  the  use  of  Brillouin  scattering  for  the  determination  of  sound  velocity  in  pure  water  were  by  Chiao,  et  al}^  and 
by  Benedek,  et  alP  Since  then,  considerable  effort  with  CW  lasers  has  been  reported,^  but  no  definitive  results  have 
previously  been  performed  using  a  pulsed  laser.  To  our  knowledge,  the  only  such  published  results  show  a  FWHM  of  the 
Brillouin  lines  of  about  3  GHz;  this  indicates  that  the  stability  and  linewidth  of  the  l^r  transmitter  were  not  sufficient^ 

Our  experimental  set  up  for  the  Brillouin  spectrum  measurements  is  shown  in  Fig.  1.  A  frequency  doubled,  injection 
seeded,  Q-switched,  pulsed  Nd:YAG  laser  provides  75  mJ  at  532  nm,  with  a  repetition  rate  of  lOHz.  Using  an  injection 
locking  technique  developed  in  our  laboratory,  it  operates  at  a  single  frequency  and  has  stable  Fourier  transform  limited 
pulses  of  10  ns,  i.e.  a  frequency  bandwidth  of  about  50  MHz.  A  horn  shaped  water  ceU  is  used  to  avoid  backreflections. 
An  additional  light  baffle  consisting  of  two  pierced  parallel  plates  was  added  in  the  bottom  of  the  water  cell  to  further 
suppress  parasitic  backreflections.  Noise  due  to  the  reflection  from  the  surface  of  the  water  is  avoid  by  using  another  thin 
baffle  plate.  The  backscattered  light  is  collected  at  an  angle  of  175^  relative  to  the  laser  beam  direction.  The  acceptance 
angle  of  the  (q>tical  system  is  *>  1*  and  is  determined  by  the  focal  length  and  the  size  of  lens  LI;  this  accq)tance  angle  does 
not  signifrcantly  increase  the  FWHM  of  the  Brillouin  lines.  The  most  critical  point  in  order  to  obtain  good  frequency 
discrimination  is  to  [vovide  a  very  weU  collimated  beam  to  the  input  of  the  Fabry-Perot  (FP)  etalon.  A  spatial  filter 
consisting  of  2  lenses  (LI  and  L2),  a  pinhole  (PI)  of  10  pm  diameter,  and  one  iris  (I)  is  used  for  this  purpose.  To  avoid 
background  radiation,  the  output  of  the  etalon  is  focused  through  pinhole  P2  onto  a  photomultiplier,  PM,  and  the  overall 
receiver  system  is  enclosed  in  a  black  box. 

The  Fabry-Perot  has  a  resolution  of  0.25  GHz  (R=%%)  and  a  Free  Spectral  Range  (FSR)  of  18  GHz.  The  etalon 
piezoelectric  transducers  are  driven  by  a  Burleigh  RC42  ramp  generator.  The  measurements  are  performed  with  a 
scanning  ramp  duration  of  2  minutes  per  FSR.  The  detector  is  a  Hamamatsu  R446  photomultiplier  tube  operating  in  the. 
analog  mode.  The  signal  is  integrated  and  av^ged  over  10  laser  shots  by  a  gated  integrator  and  boxcar  averager.  It  is 
digitized  and  stored  using  a  Labview  program  on  a  MacII  computer. 

Fig.  2  shows  a  measured  frequency  spectrum,  superposed  on  the  theoretical  one  for  a  water  temperature  of  29‘C 
(corresponding  to  a  Brillouin  shift  of  7.55  GHz).  Only  one  free  spectral  range  of  the  Fabry-Perot  is  presented  in  the  figure. 
The  experimental  points  are  raw  data  from  a  single  scan;  no  further  iHxx:essing  has  been  p^cumed.  There  is  very  good 


agreement  with  thewy.  The  strong  central  peak  is  due  to  particulate  contamination  of  the  water,  which  is  apparently  due 
mostly  to  oxidation  of  the  metal  light  baffle  submerged  in  the  horn  shaped  water  cell. 


Fig.  1.  Brillouin  spectrum  measurement  set  up.  LI,  L2,  and  L3  are  biconvex  lenses  of  focal  length  50  cm,  20  cm  and  8 
cm,  PI  is  a  10  |im  diameter  pinhole,  I  is  an  iris  used  to  select  the  lower  interference  order,  P2  is  a  1  mm  diameter 
pinhole,  SHG  is  second  harmonic  generator,  M  are  mirrors,  and  PM  is  the  photomultiplier. 
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Fig.  2.  Experimental  (dots)  and  theoretical  Brillouin  spectrum  of  water  at  29*C  using  the  second  harmonic 
of  a  pulsed  injection  seeded  Nd;  YAG  laser  (532  nm). 


3.  BRILLOUIN  SHIFT  MEASUREMENT 


As  discussed  in  the  introduction,  a  Brillouin  LIDAR  based  on  a  version  of  the  edge  technique  is  under  development. 
Using  the  edge  technique,  atmospheric  measurements  at  1064  nm  have  shown  wind  velocity  accuracies  of  0.2S  m/s, 
corresponding  to  frequency  changes  smaller  than  470  kHz.^^  This  is  the  most  sensitive  method  yet  demonstrated  for  small 
frequency  shift  LIDAR  measurements.^^  V^ions  of  the  edge  technique  had  already  been  considered  by  Hirschberg,  et  al? 
and  by  Hickman,  et  al}  for  oceanic  Brillouin  shifts  measurements  using  a  laser  wavelength  of  S32  nm.  In  the  latter  study 
the  statistical  uncertainty  8vb  in  the  measurement  of  the  Brillouin  shift  is  calculated  assuming  a  single  laso'  shot  with  an 
energy  of  500  mJ  and  a  depth  resolution  of  1.1  meter.  For  y  <  0.04  m'^  shot  noise  limited  uncertainties  in  5vb  of  less  than 
1  MHz  can  be  expected  for  depths  down  to  75  meters,  and  of  less  than  4  MHz  down  to  100  metos.  These  correspond  to 
sound  velocity  measurement  accuracies  that  can  be  better  than  0.3  m/s  and  1  m/s}  respectively. 

We  are  developing  a  version  of  the  edge  technique  that  makes  use  of  the  edges  of  absorption  lines  in  the  molecular 
spectra  of  I2  and/or  Br2,  as  reposed  to  the  Fabry-Perot  interferometers  used  in  other  applications'"^’  ®  of  the  edge  technique. 
The  schematic  of  such  a  LIDAR  receiver  is  shown  in  Fig.  3.  The  second  harmonic  of  an  injection-seeded  pulsed  Nd:YAG 
lasCT  is  used  as  the  transmitter.  The  lasCT  signal  backscattered  by  the  ocean  water  is  collected  by  a  telescope,  and  is  first 
directed  through  an  atomic  cell  (used  as  a  blocking-frlter)  where  the  central  line  of  the  backscattered  frequency  spectrum  is 
absorbed.  The  backscattered  signal  transmitted  by  the  blocking  filter  is  then  split  into  two  equal  parts;  one  part  is  detected 
to  produce  a  normalization  signal  (S2).  The  other  part  passes  through  the  edge  filter  and  is  detected  to  produce  the  signal 
(SI).  The  Nd:YAG  laser  frequency  stabilization  is  performed  by  picking  off  a  small  part  of  the  cw  seed  laser  beam, 
frequency  doubling  it  in  an  external  ring  build-up  cavity,  and  locking  the  seed  laser  to  the  absorption  line  of  the  blocking- 
filter.  The  return  signals  SI  and  S2  are  sampled  every  2-3  ns  to  produce  a  range-resolved  measurement  (of  course,  the 
temporal  resolution  is  limited  by  the  10  nanosecond  laser  pulse). 
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Fig.  3.  Brillouin  LIDAR  receiver  schematic.  The  Brillouin  shift  is  determined  from  the  ratio  S1/S2.  BS  are  beam  splitters,  M  are  mirrors, 
SHG  is  second  harmonic  generators,  M  are  mirrors,  and  PD  are  photo  detectors.  Solid  lines  correspond  to  light  paths,  and  dotted 
lines  to  electrical  connections. 

As  an  example  of  an  edge  filter,  consider  a  '^12  absorption  cell.  Fw  a  laser  wavelragth  of  532.573806  nm,  it  provides 
two  apprc^riately  positioned  absorption  edges  (see  Fig.  4).  Because  of  the  steep  average  slope  of  the  absorption  edges, 
small  changes  in  the  Brillouin  shift  produce  large  changes  in  the  transmission  of  the  edge  filter.  The  temperature  of  the  gas 
and/or  the  vapor  pressure  in  the  cells  can  be  changed  for  a  fine  adjustment  of  the  absorption  line  slopes  and  frequencies. 

An  important  advantage  of  this  edge  filler  version  is  that  a  change  in  the  incident  laser  frequency  ot  any  Doppler  shifts 
due  to  movements  will  increase  the  edge  filter  signal  OMresponding  to  one  Brillouin  line  and  decrease  Ae  signal 
corresponding  to  the  other;  thus  such  systematic  effects  are  canceled  to  first  order.  The  ratio  S1/S2  provides  a 


normalization  which  makes  the  frequency  discrimination  independent  of  the  intensity  of  the  light  backscattered  by  the 
ocean  water,  i.e.  independent  of  the  attenuation  due  to  water  absorption  or  hydiosols  scatt^ng. 

An  ideal  blocking  filter  should  exhibit  total  absorption  for  the  central  line  of  the  backscattered  spectrum,  and  good 
transmission  for  the  Brillouin  peaks.  The  FWHM  of  the  central  line  is  mainly  the  linewidth  of  the  laser  transmitter,  and  is  « 
50  MHz.  Fig.  5  shows  the  absorption  spectrum  of  Br2  superposed  on  theoretical  Brillouin  lines  with  a  shift  of  7.5  GHz, 
and  for  the  same  incident  laser  wavelength  as  for  the  I2  data  in  Fig.  4.  This  absorption  line  provides  a  satisfactory, 
although  not  perfect  blocking  filter.  Of  course,  the  convolution  of  the  Brillouin  lines  and  of  the  transmission  spectrum  of 
Br2  has  to  be  c(xisidered  as  the  input  to  the  edge  frlter. 


Fig.  4.  Edge  filter:  theoretical  Brillouin  lines  for  a  frequency 
shift  of  7.5  GHz  superposed  on  the  absorption  spectrum 
of  Small  changes  of  the  Brillouin  shift  make  large 
changes  in  the  transmitted  signal.  The  incident 
wavelength  for  the  Brillouin  spectrum  is  532.573806  nm. 


Fig.  5.  Blocking  filter  :  The  ideal  blocking  filter  should  have  a 
good  transmission  for  the  Brillouin  lines  and  a  total 
absorption  for  the  central  line  of  the  backscattered 
spectrum.  The  absorption  lines  of  ^^Br2.  provide  a 
satisfactory  blocking  filter  for  an  incident  laser 
wavelength  of  532.573806  nm. 


4.  PROSPECTS  AND  CONCLUSION 


We  have  shown  that  an  injection  seeded  pulsed  laser  can  be  used  to  obtain  completely  resolved  measurements  of 
Brillouin  lines.  This  verifies  that  Brillouin  shift  data  can  be  used  to  measure  sound  velocity  as  a  function  of  depth  in  the 
ocean.  A  version  of  the  edge  technique  for  the  measurement  of  the  Brillouin  shift  has  been  outlined. 

In  comparison  to  other  LIDAR  receiver  designs  (using,  for  example,  a  Fabry-Perot)  the  I2  and/or  Br2  absorption  cells 
provide  a  simpler  and  m«e  robust  tqiivoach  that  requires  no  critical  alignments.  The  optical  efficiency  of  a  receiver  using 
the  abscMption  lines  discussed  here,  would  be  as  good  as  that  using  a  Fabry-Perot  based  receiver.  A  systematic  search  for 
other  combinations  of  transmitter  wavelength  and  absorption  lines  is  underway.  This  includes  direct  experimental 
examination  and  comparison  (rf  the  absorption  lines  (rf  ^27i2  l29i2^  793^2,  and  *^Br2. 
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Abstract 

In  this  paper,  we  rigorously  analyze  Brillouin  LIDAR  measurements 
of  the  temperature  and  sound  velocity.  There  are  five  variables: 
Brillouin  shift,  sound  velocity,  index  of  refraction,  temperature  and 
salinity  and  three  well-known  relations  between  them.  In  order  to 
express  all  of  them  in  terms  of  a  single  independent  variable  (e.g. 
Brillouin  shift),  one  more  relation  is  required.  For  the  purpose  of  such 
a  relation,  the  present  study  employs  the  mean  value  and  standard 
deviation  of  salinity  data  obtained  by  analyzing  a  very  large  set  of 
historical  in-situ  data.  Temperature  and  sound  velocity 
measurements  can  be  expected  down  to  100  meters  depth  to  have  an 
accuracy,  as  high  as  of  OJ  degree  and  0.2  m/s,  respectively. 
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1.  Introduction 


Accurate  remote  sensing  of  the  temperature  and  soimd  velocity 
structure  of  the  upper-ocean  mixed  layers  is  of  major  importance  to 
the  understanding  of  the  physicai  and  biologicai  behavior  of  the 
ocean.1'  2  There  are  also  many  other  applications;  it  is,  for  example,  a 
necessary  input  to  weather  forecast  and  climate  change  studies.^ 
There  currendy  exists  no  instrumentation  for  remote  sensing  of  these 
parameters,  and  only  in’SJtu  measurements  have  been  performed 
with  the  required  accuracy.  However,  in-situ  approaches  do  not  allow 
the  rapid,  accurate,  and  real-time  range-resolved  monitoring  that  is 
needed.  UDAR  remote  sensing  provides  a  promising  solution.  Raman 
and  Brillouin  backscattering  are  both  sensitive  to  temperature,  and 
Brillouin  scattering  is  directly  related  to  sound  velocity.  Both 
processes  have  been  previously  considered  for  UDAR  measurements 
of  the  vertical  structure  of  temperature  and  sound  velocity  in  the 
ocean.2*i3  Major  relative  limitations  of  the  Raman  method  inciude:  (1) 
The  large  Raman  frequency  shift,  which  means  the  transmitted  laser 
field  and/or  the  backscattered  field  must  be  in  a  relatively  highly 
absorbing  region  of  the  water  spectrum.  (2)  The  Raman 
backscattering  spectrum  has  a  much  wider  frequency  distribution 
than  the  Brillouin  backscattering  one.  Consequently,  the  Raman 
receiver  must  have  a  much  wider  bandwidth  and  is  therefore  much 
more  susceptible  to  background  signals  and  noise.  As  a  result  of 
these  limitations,  the  Raman  technique  does  not  appear  promising. 

Several  studies  have  shown  the  potential  accuracy  of  UDAR 
measurements  of  sound  velocity  and  temperature  in  the  ocean  by 
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measuring  the  spectral  shift  of  the  backscattered  Brillouin  lines.2-6.8 
Leonard  and  Sweeney^^  have  even  said  that  to  achieve  a 
temperature  accuracy  in  the  range  of  millidegrees  over  two  or  three 
attenuation  lengths,  advanced  Brillouin  methods  must  be  invoked. 
However,  taking  into  accovmt  other  relevant  factors,  we  show  in  this 
paper  that  for  Brillouin  measurements  alone  (i.e.  no  additional 
complementary  measurements),  the  expected  accuracy  is  only  some 
tenths  of  a  degree;  this  is  however  sufficient  for  many  applications. 
Previously,  technical  limitations  had  prevented  the  realization  of 
such  Brillouin  LIDAR  instrumentation.  Nevertheless  the  required 
technology,  in  particular  pulsed  lasers  with  a  Fourier  transform 
limited  bandwidth,  are  now  commercially  available,^  and  a  receiver 
based  on  an  edge  technique^'^  using  absorption  lines  in  Brz  and  I2  now 
appears  feasible.^^  We  emphasize  that  the  objective  of  this  paper  is 
to  determine  general  limits  on  the  accuracy  of  UDAR  Brillouin  shift 
measurements  and  not  to  describe  an  overall  system. 

The  Brillouin  lines  consist  of  two  inelastically  scattered 
Lorentzians  centered  symmetrically  at  ±vb  with  respect  to  the 
transmitted  laser  line;  for  water  the  FWHM  of  the  Brillouin  lines  is 
approximately  0.5  GHz  and  vg  is  7  to  8  GHz.  For  a  given  incident  laser 
wavelength  A,  the  relation  between  the  Brillouin  shift  vg  and  the 
velocity  of  sound  vs  in  ocean  water  is  given  for  180*  backscattering 
by,^-^ 

VB(S,T.p)  =  ^^^^V3(aT.p),  (1) 


On  Brillouin  UDAR ... . 


3/36 


Fry,  et.  al. 


where  n  is  the  refractive  index  of  water,  S  is  the  salinity,  T  is  the 
temperature,  and  p  is  the  pressure. 

To  derive  the  sound  velocity  from  the  Brillouin  shift 
measurement,  the  refractive  index  n  is  has  frequently  been  assumed 
to  be  known.^'  *•  A  more  accurate  approach  had  been  considered  by 
ffickman,  et  al,y^  who  point  out  that  the  uncertainty  in  our  knowledge 
of  the  refractive  index  due  to  our  uncertainty  in  T  and  S  contributes 
to  the  error  in  the  measurement  of  the  sound  velocity.  Using  an 
empirical  relation  for  the  refractive  index  as  a  function  of  T  and  $,^6 
and  assuming  that  T  is  known  to  ±2  *C,  and  S  is  known  to  ±l%o,  they 
estimate  that  the  refractive  index  variability  limits  the  sound 
velocity  accuracy  to  approximately  0.5  m/s. 

To  derive  the  temperature  from  the  Brillouin  shift  measurement, 
an  empirical  equation  is  generally  used  to  describe  the  dependence 
of  the  sound  velocity  on  S,  T,  and  p.^^  Then  using  the  measured 
Brillouin  shift  and  assuming  knowledge  of  the  salinity,  the  refractive 
index,  and  the  pressure,  Eq.  1  gives  the  temperature.^* 

In  this  paper,  we  will  provide  a  rigorous  approach  to  these 
interdependencies.  We  will  identify  five  variables  (vg,  S,  T,  vs,  n)  and 
since  we  will  have  three  relations  between  them,  two  of  these 
variables  must  be  independent,  and  three  must  be  dependent.  Since 
we  will  measure  vg,  we  take  it  as  one  of  the  independent  variables. 
We  take  S  as  the  other  independent  variable,  and  obtain  empirical 
equations  for  the  dependent  variables  T,  vs,  and  n  as  a  function  of 
the  independent  variables  vg  and  S.  Such  equations  for  T  and  vs  as  a 
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function  of  vq  and  S  will  be  given  in  this  paper.  Throughout  the 
discussion  we  will  neglect  the  dependence  on  pressure  since  it  is  a 
well-known  function  of  depth  and  can  be  included  later. 

Since  S  is  the  other  independent  variable,  a  very  large  set  of  in- 
situ  salinity  meastirements  was  analyzed  to  determine  the  salinity 
mean  value  and  standard  deviation  as  a  function  of  depth  and 
location.  These  results  enable  us  to  determine  a  limit  on  the  accuracy 
in  the  retrieval  of  T  and  vs  from  measurements  of  the  Brillouin  shift, 
when  historical  data  are  used  for  the  salinity  mean  value  and 
standard  deviation. 

The  relations  between  vg,  S,  T,  vs  and  n  are  examined  in  Sect.  2. 
The  limit  on  the  accuracy  of  measurements  of  T  and  vs  are  discussed 
in  Sect.  3;  effects  due  to  the  estimates  of  the  variability  in  S  and  to 
the  errors  in  the  measurement  of  vg  are  included. 

2.  Relations  between  Brillouin  shift,  salinity,  temperature, 
and  sound  velocity  in  the  ocean. 

The  most  widely  used  empirical  equation  for  the  dependence  of 
sound  velocity  on  salinity,  temperature,  and  pressure  in  natural 
water  is  given  by:^^ 

Vs(S,T,p)  =  Co+CiT+C2T2+C3T3+C4S+C5S^+cJS+C7T^S+f(S,T,p),  (2) 

where  the  coefficients  ci  are  constants  and  f(S,T,p)  represents  all  the 
terms  containing  the  presstire.  The  ci's  with  T  in  ’C,  S  in  %o,  p  in  atm. 
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and  Vs  in  m/sec  are  given  in  Table  L  Typical  values  for  the  present 
application  range  between  0  “C  and  30  “C  for  T,  0  %o  (fresh  water)  and 
40  %o  for  S,  and  1  atm  and  10  atm  for  p  (corresponding  to  depths  of  0 
and  approximately  100  meters,  respectively).  Typical  values  for  the 
changes  in  vs  due  to  changes  in  p,  T,  and  S  are: 


'^]  -0.162^, 

^  ^  Jsmsn,,  T»io*c,  p-o 


„.m/sec 
=  3.5  -,and 

S-36V  T»WC.  p-0  ^ 

=1.22=^. 

V  Jsm36%,,  T*wc.  p=0  ^ 


(3) 


It  should  also  be  noted  that  the  horizontal  variability  of  pressure 
at  a  given  depth  is  negligible.  For  simplicity,  the  pressure  is  assumed 
to  be  zero  for  the  present  calculations,  i.e.  f(S,T,p)=0.  Generalization  to 
other  pressures  is  straightforward. 


An  empirical  equation  for  the  dependence  of  the  refractive  index 
n  of  sea  water  at  atmospheric  pressure  on  temperature,  salimty,  and 
wavelength  is  given  by:^* 


where  Xis  the  wavelength  in  nanometers,  T  is  in  *C,  S  is  in  %o,  and  the 
coefficients  nj  are  constants,  given  in  Table  H. 
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Again,  the  five  variables  of  interest  are  vg,  S,  T,  vs,  and  n.  Since 
Eqs.  1, 2,  and  4  provide  three  relations  between  them,  any  two  of  the 
five  can  be  chosen  as  independent  Our  choices  are  vg  and  S  since  the 
former  is  our  measurement  result,  and  the  latter  is  the  most  stable 
and  well-known,  as  we  will  see  in  section  3.1.  With  these  two  as 
independent  variables,  the  other  three  variables  can  be  expressed  as 
T(S,vg),  vs(S,vg),  and  n(S,vg). 

To  obtain  empirical  equations  for  T(S,vg)  and  vs(S,vg),  we  use 
Eqs.  1,  2,  and  4  to  generate  a  table  containing  vs(S,T),  n(S,T),  and 
vg(S,T)  at  fixed  values  of  S  and  T  (Table  III).  Then,  taking  the  entries 
for  S,  vg,  and  T,  a  least-squares  fit  is  performed  to  determine  a 
function  T(S,vg)  of  the  form: 

T(S,  Vb)  =  to  +  t^CVe  -  7.5) + tj(vB  -  7.5)*  +  %(yj,  -  7.5)®  +  t^Cvg  -  7.5)*  + 

S[tB  +  te(vB  -  7.5)  +  t^CvB  -  7.5)*  +  tgCvB  -  7.5)®]. 

The  resulting  tj  are  given  in  Table  IV.  Numerical  evaluations  of 
Eq.  5  are  shown  in  Fig.  1,  where  T  is  plotted  as  a  ftmction  of  S  for 
various  values  of  vg  (7.4  GHz  to  7.7  GHz  in  steps  of  50  MHz),  and  in 
Fig.  2,  where  T  is  plotted  as  a  function  of  vg  for  various  values  of  S 
(30%o  to  40%o  in  steps  of  2.5%o). 

Similarly,  for  the  sound  velocity,  we  use  the  entries  for  S,  vg,  and 
Vs,  to  perform  a  least-squares  fit  and  obtain  a  function  vs(S,vg)  of  the 
form: 

Vs(S,Vb)  =  c',+c',  (Vb  -7.5)+c',(Vs  -7.5)“  +c'a(v,-7.5)“  + 

S[c'.  +c'.  (V,  -  7.5) + c',  (V,  -  7.5)“} 
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The  resulting  c'j  are  given  in  Table  V.  Numerical  evaluations  of 
Eq.  6  are  shown  in  Fig.  3,  where  vs  is  plotted  as  a  function  of  S  for  vq 
from  7.4  to  7.7  GHz  in  steps  of  50  MHz,  and  in  Fig.  4,  where  vs  is 
plotted  as  a  function  of  for  S  from  30%o  to  40%o  in  steps  of  5%o.  Eqs. 
5  and  6  together  with  the  coefficients  in  Tables  m  and  IV  are  valid 
for  30  %o  <  S  <  40  %o  and  7.0  GHz  <  vg  <  8.0  GHz. 

Fig.  1  shows  that  the  temperature  has  a  relatively  weak 
dependence  on  the  salinity.  Specifically,  from  Eq.  5,  we  find 

=  t5  =  -0.4  ®C/*o.  (7) 

v,«7.5GHs 

For  example,  if  the  uncertainty  in  the  salinity  is  about  0.5%o,  the 
uncertainty  in  the  temperature  measurement  is  ^0.2  *C.  We  will  see 
in  section  3.1  that  such  an  uncertainty  in  the  salinity  of  this 
magnitude  or  less  occurs  frequendy  in  historical  data  for  salinity  in 
the  ocean  at  a  given  depth  and  location. 

On  the  other  hand.  Fig.  2  shows  a  relatively  strong  dependence  of 
T  on  vg.  From  Eq.  5  we  have 

=  + St,  =  0.056  «C/MHz.  (8) 

S-36«>,vi-7^Ki 

Uncertainties  in  Brillouin  shift  measurements  of  1  MHz  would  lead  to 
temperature  uncertainties  of  about  0.06  *C.  As  will  be  discussed  in 
section  3.2,  such  a  measurement  accuracy  can  be  expected  with  a 
Brillouin  IIDAR. 


3T 

OVg 


dS 
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The  weak  dependence  on  the  salinity  is  even  more  pronounced 
in  the  data  of  Fig.  3  for  the  sound  velocity  as  a  function  of  the 
salinity.  Specifically,  Eq.  6  gives 


Va-7^Hi 


=  -0.24  Tnlsl%c; 


thus  an  uncertainty  in  the  salinity  of  0.5%o  gives  an  uncertainty  in 
the  sound  velocity  of  «  0.12  m/s.  Actually,  a  very  rough  knowledge 
of  the  salinity  (a  few  %o)  is  sufficient  to  determine  the  sound  velocity 
with  an  imcertainty  less  than  1  m/s. 


Fig.  4  shows  that  the  dependence  of  vs  on  vq  is  relatively  strong. 
Eq.  6  gives: 


^Vs 

3vb 


=  c'l  +Sc'5  =  0.2  m/8 / MHz. 


(10) 


IS»36«*,Va»7.5GHs 


Thus,  a  measurement  error  in  the  Brillouin  shift  of  1  MHz 
corresponds  to  an  uncertainty  of  0.2  m/s  for  the  sound  velocity. 


Compared  to  previous  studies,  use  of  the  three  relations  (Eqs.  1, 
2,  and  4)  has  enabled  us  to  rephrase  the  problem  so  that  only  the 
uncertainties  in  vg  and  S  contribute;  for  example,  temperature  has 
become  a  derived  quantity  and  its  uncertainty  is  irrelevant  to  the 

determination  of  the  uncertainty  in  Vs. 
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3.  Measurement  uncertainty  in  the  temperature  and  sound 
velocity. 

The  root  mean  square  error  (uncertainty)  in  measurements  of 
the  temperature,  5T,  and  sound  velocity,  5vs,  can  be  expressed  in 
terms  of  the  uncertainties  in  the  salinity  S  and  Brillouin  shift  5vb: 


The  following  two  sections  discuss  estimates  for  the  uncertainty 
in  the  salinity  and  in  the  Brillouin  shift  measurement 

3.1  Variability  of  the  salinity  in  the  ocean. 

A  standard  deviation  5S  for  the  salinity  has  been  estimated  by 
referring  to  the  Oceanographic  Station  Profile  Time  Series.^^  It 
provides  a  compilation  of  oceanographic  data  obtained  from 
repetitive  samples  taken  along  various  ocean  sections  and  at  fixed 
coordinates  (oceanic  weather  stations)  between  1900  and  1990.  We 
have  analyzed  this  data  set  for  the  salinity  variability  at  113 
geographical  locations,  distributed  over  all  the  oceans  of  the  earth 
(more  than  75,000  measurements).  Each  location  corresponds  to  an 
area  of  about  500x500  km^.  For  each  location  we  have  computed  the 
mean  value  and  the  standard  deviation  at  various  depths.  A  typical 
example  is  shown  in  Fig.  5  at  32“15'  S  17*30'  E,  near  the  Cape  of  Good 
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Hope,  at  a  depth  of  10  meters,  and  in  an  area  of  30x100  km2.  The 
variability  of  the  salinity  is  constant  over  a  period  of  25  years. 

The  results  of  our  analysis  are  summarized  in  Fig.  6,  where  the 
cumulative  distribution  of  the  salinity  standard  deviation  is  shown 
for  depths  of  10  meters  (Fig  6.a)  and  100  meters  (Fig  6.b).  It  is 
expressed  as  a  fraction  of  the  total  number  of  locations.  In 
particular,  for  a  depth  of  10  meters,  50%  of  the  locations  have  a 
salinity  standard  deviation  5S  less  than  0.5%o,  and  in  92%  of  the 
cases,  5S  is  less  than  l%o.  For  a  depth  of  1(X)  meters,  8S  is  less  than 
0.2  %o  for  43%  of  the  locations,  and  is  less  than  l%o  for  99%  of  the 
locations.  As  expected,  the  standard  deviation  generally  decreases 
with  depth. 

Some  locations  exhibit  seasonal  variations,  which  could  be  taken 
into  accoimt  in  the  reduction  of  8S.  An  example  is  shown  in  Fig.  7  for 
measurements  performed  in  the  middle  of  the  North  Pacific  (34*00'  N 
164*00'  E)  at  depths  of  10  m,  50  m,  and  100  m,  in  an  area  of  50x50 
km2.  For  this  location,  8S  would  be  significantly  smaller  if  the 
seasonal  trend  is  considered. 

We  have  also  noticed  that  for  a  given  station,  8S  does  not 
substantially  change  when  increasing  or  decreasing  the  horizontal 
sampling  area  or  when  increasing  the  thickness  of  the  water  cell 
considered  at  a  given  depth.  This  shows  that  the  salinity  is  fairly 
stable,  and  that  historical  data  provide  very  useful  salinity  data  that 
can  be  used  to  determine  T  and  vs  from  Brillouin  shift 

measurements. 
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3.2.  Brillouin  shift  measurement  accuracy. 

The  use  of  injection  locking  techniques  enable  one  to  produce 
stable  laser  pulses  having  a  Fourier  transform  limited  bandwidth  of 
approximately  50  MHz  at  a  wavelength  of  532  nm.20'  21  jjie 
linewidth  of  the  pulsed  laser,  as  well  as  its  stability,  are  no  longer  a 
technical  limitation  for  measurements  of  the  Brillouin  shift. 

For  a  range-resolved  measurement,  a  non-scanning  high- 
resolution  spectroscopic  technique  is  needed.  The  so-called  edge 
technique^'^  has  been  demonstrated  to  be  a  very  powerful  method  for 
the  measurement  of  small  frequency  changes  in  an  atmospheric 
Doppler  LIDAR  application.  Laboratory  measurements^^  have  given 
wind  velocity  with  an  uncertainty  of  0.25  m/s,  corresponding  to  a 
frequency  change  of  less  than  470  kHz  at  1064  nm.  This  is  presently 
the  most  sensitive  method  to  be  demonstrated  for  small  frequency 
shift  UDAR  measurements.23  Such  a  technique  has  been  considered 
for  oceanic  Brillouin  shift  measurements  using  a  laser  wavelength  of 
532  nm.2’ 

In  a  recent  study^  the  statistical  uncertainty  Svg  in  the 
measurement  of  the  Brillouin  shift  is  calculated  for  various  effective 
attenuation  coefficients  7  of  water,  assuming  a  single  laser  shot  (500 
mj),  a  depth  resolution  of  1.1  meter,  and  known  technical 
parameters.  For  y  <  0.04  m-i  one  can  expect  5vb  to  be  less  than  4  MHz 
for  depths  down  to  100  meters,  and  less  than  1  MHz  for  depths  down 
to  75  meters.  Based  on  these  results,  we  consider  two  typical 


Brillouin  shift  measurement  errors  Svg  of  1  MHz  and  4  MHz  in  the 
following  estimates. 

3.3.  Temperature  and  sound  velocity  measurement 
accuracy. 

Eqs.  5  and  6  and  the  estimates  of  aS  and  Svg  made  in  the  last  two 
sections  allow  the  evaluation  of  the  temperature  and  sound  velocity 
measurement  errors  given  by  Eqs.  11  and  12. 

Fig.  8  shows  the  standard  deviation  in  the  temperature  5T  for 
5S=0.5%o  and  Svg  =  1  MHz.  This  uncertainty  8T  is  plotted  as  a  function 
of  S  for  typical  oceanic  values  between  30%o  and  40%o,  and  with  vg  as 
a  parameter  from  7.4  GHz  to  7.7  GHz  in  steps  of  50  MHz.  For  a  typical 
value  of  vb=7.5  GHz,  the  temperature  resolution  is  about  0.2  “C. 

Fig.  9  shows  5T  (left  side  y-axis)  as  a  function  of  5S  superposed 
on  the  cumulative  distribution  for  the  salinity  standard  deviation  at  a 

depth  of  100  meters  (right-side  y-axis)  as  in  Fig.  6b.  For  these  results 
we  have  taken  vg=7.5  GHz  and  S=35%o.  When  5S  is  smaller  than 

«0.3%o,  the  uncertainty  Svg  is  the  dominant  contribution  to  5T.  If  8vb 
is  less  1  MHz,  the  uncertainty  in  the  temperature  is  less  than 
0.15  *C  for  more  than  50  %  of  the  locations.  At  99%  of  the  locations 
5S<l%o,  thus  if  6vb-4  MHz,  5T  is  less  than  0.5  *C  at  these  locations. 

Fig.  10  shows  the  standard  deviation  of  the  soimd  velocity  Sv*  for 
5S=0.5%o  and  uncertainty  5Vs  is  plotted  as  a  function 

of  S  for  30%o  ^  S  ^  40%o  and  7.4  GHz  ^  vg  ^  7.7  GHz  in  steps  of  50 
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MHz.  Notice  that  throughout  the  range  of  parameters  for  Fig.  10,  the 
sound  velocity  imcertainty  is  always  less  than  about  0.25  m/s. 

Fig.  11  shows  8vs  (left  side  y-axis)  as  a  function  of  8S  superposed 
on  the  cumulative  distribution  of  the  salinity  standard  deviation 
(right  side  y-axis)  at  a  depth  of  100  meters,  as  in  Fig.  6b.  For  these 
results  we  have  also  taken  vg*7.5  GHz  and  S“3596o.  The  small  slope  of 
Svs(SS)  is  a  consequence  of  the  weak  dependence  of  vs  on  S.  The 
Brillouin  frequency  shift  is  the  dominant  contribution  to  the 
imcertainty  in  a  sound  velocity  measurement.  Assuming  5vb  is  less 
than  1  MHz,  8vs  is  less  than  0.3  m/s  for  more  than  99%  of  the 
locations.  Even  for  8vb=4  MHZ,8v,  is  still  always  less  than  1  m/s. 

In  fresh  water  there  are  only  4  parameters  since  S  is  already 
known  to  be  zero.  Since  there  are  three  relations,  measurement  of 
one  parameter,  v^,  is  sufficient  to  determine  the  others.  For  the 
determination  of  the  accuracy  of  Brillouin  LIDAR  measurements  in 
fresh  water,  we  set  S  and  8S  to  0.  Then  for  8vb=1  MHz,  we  find  8T  is 
about  0.06  "C,  and  8vs  is  about  0.2  m/s. 

Finally,  it  should  be  noted  that  if  one  more  of  the  variables  vg,  S, 
T,  vs»  n,  or  equivalently  one  more  relation  between  them  is  also 
measured,  then  all  the  others,  including  salinity  can  be  determined 
without  the  necessity  of  recourse  to  historical  data  for  salinity.  There 
are  several  possibilities,  including  simultaneous  use  of  Raman 
scattering,  or,  for  example,  measurements  of  conductivity.  For  some 
situations  there  is  also  an  S-T  relationship”  that  may  be  applicable  to 
increase  the  measurement  accuracy. 
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4.  Conclusions 


We  have  obtained  rigorous  relations  between  vb,  S,  T,  vs,  and  n, 

and  have  used  them  to  determine  the  limits  of  accuracy  of  a  Brillouin 
UDAR  system.  We  have  also  analyzed  a  large  set  of  historical  data  to 
determine  the  mean  and  the  standard  deviation  of  the  salinity.  Use  of 
this  additional  parameter  allows  us,  to  show  that  the  temperature 
and  sound  velocity  can  be  retrieved  with  an  accuracy  sufficient  for 
many  applications.  Implementation  of  a  Brillouin  UDAR  system  is 
underway,  and  the  first  frequency  resolved  measurement  of  the 
Brillouin  lines  using  a  pulsed  laser  have  been  recently  obtained.^^ 
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Figure  Captions 


Fig.  1.  Temperature  T(\)b»  S)  shown  as  a  function  of 
salinity  S  with  the  Brillouin  frequency  shift  \)b  as  a 
parameter  between  7.4  and  7.5  GHz  in  steps  of  50 
MHz. 

Fig.  2.  Temperature  T(vb»  S)  shown  as  a  function  of  the 
Brillouin  frequency  shift  db  for  values  of  salinity  S 
from  30  %o  to  40  %o  in  steps  of  2.5  %o. 

Fig.  3.  Sound  velocity,  vs(\)b,  S),  shown  as  a  function  of 
salinity  S  with  the  Brillouin  frequency  shift  db  ^  ^ 
parameter  (in  steps  of  50  MHz). 

Fig.  4.  Sound  velocity,  vs(vb,  S),  shown  as  a  fimction  of 
the  Brillouin  frequency  shift  ub  for  values  of 
salinity  S  from  30  %o  to  40  %o  in  steps  of  5  %o. 

Fig.  5.  Variability  of  the  salinity  at  32“15'  S  17*30'  E,  in 
the  South  Atlantic,  near  the  Cape  of  Good  Hope,  at 
a  depth  of  10  meters.  Sampling  area:  30>d00  km^. 

Fig.  6.  Cumulative  histogram  of  the  variability  of  the 
salinity  standard  deviation  for  113  locations  for  a 
depth  of  10  meters  (6a)  and  100  meters  (6b). 
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Fig.  7.  Seasonal  variability  of  the  salinity  at  34*00’  N 
164*00'  E,  in  the  middle  of  the  North  Pacific,  at  a 
depths  of  10  m,  50  m,  and  100  m.  Sampling  area: 
50x50  km2.  The  standard  deviation  represents  the 
overall  variability,  the  seasonal  trend  is  not  taken 
into  accoimt  in  the  calculation. 

Fig.  8.  The  uncertainty,  8T,  in  the  determination  of 
temperature  as  a  function  of  the  salinity  for 
8S»0.5%o  and  5ub=1  MHz,  and  various  Brillouin 
shifts  from  7.4  to  7.7  GHz  in  steps  of  50  MHz. 

Fig.  9.  On  the  left-side  y-axis,  the  uncertainty,  8T,  in  the 
determination  of  temperature  is  plotted  as  a 
function  of  the  salinity  standard  deviation,  for 
8vb*1  MHz  (solid  lines)  and  5vb=4  MHz  (dotted 
lines).  In  both  cases  we  have  taken  ub=7.5  GHz  and 
S=35  %o.  It  is  superposed  on  the  salinity  standard 
deviation  cumulative  probability  distribution  at 
100  meters,  on  the  right-side  y-axis  (same  as  Fig. 
6b). 

Fig.  10.  The  imcertainty  5vs  in  the  determination  of  sound 
speed  as  a  function  of  salinity  with  ub  as  a 
parameter  that  is  varied  from  7 .4  GHz  to  7 .7  GHz  in 
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steps  of  50  MHz.  For  these  calculations  we  have 
taken  5ub“1  MHz  and  $=0.5%o. 

Fig.  11  On  the  left-side  y-axis,  the  uncertainty,  5vs,  in  the 
determination  of  sound  velocity  is  plotted  as  a 
function  of  the  salinity  standard  deviation,  for 
5^8==  1  MHz  (solid  lines)  and  5t)B=4  MHz  (dotted 
lines).  In  both  cases  we  have  taken  t)B=7.5  GIfe  and 
S=35  %o.  It  is  superposed  on  the  salinity  standard 
deviation  cumulative  probability  distribution  at 
100  meters,  on  the  right-side  y-axis  (same  as  Fig. 
6b). 
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Tables 


Table  I.  Coefficients  in  the  empirical  expression  for 
the  sound  velocity  vs(S,T).^^ 


Co»1402.392 

ci=5.01 109398873, 

C2=-0.0550946843172 

C3=0.00022153596924, 

C4=1.32952290781 

C5=0.0001289557568, 

C6=-0.0127562783426 

C7=0.0000968403  15641. 

Table  II.  Coefficients  in  the  empirical  expression  for  the  refractive  index 

n(S.T,X).'» 


no=1.31405, 

ni=1.779  X  10-4 

n2=-1.05  X  10-^, 

n3*1.6  X  10‘*, 

n4=-2.02  X  10-6, 

n5=15.868. 

n6=0.01155, 

n7=-0.00423, 

n8*-4382. 

nq-1.1455  X  10^: 
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Table  III.  n,  vs  [m/s],  vb  [GHz],  for  a  range  of  values  of  30%o  <  S  <  40%o  and  0“C  <  T  <  40"C; 
p>-0  atm,  and  ^532.57  nm. 
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Table  IV.  Coefficients  in  the  empirical  expression  for  T(S,  vb) 

at  532.57  nm. 


to=23.689 

ti=68.696 

t2=5 1.939 

t3»225.809 

U-274.492 

t5=-0.40554 

t6=-0.35856 

t7="0.56583 

t8=-5.4252 

Table  V.  Coefficients  in  the  empirical  expression  for  vsIS^Vb) 

at  532.57  nm. 


c'o=1496.17 

c'i=205.08 

c'2=30.924 

c’3=35.242 

c'4=-0.2406 

c'6="0.6336 

c'5=-0.08531 
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Fig.  1.  Temperature  T(ub,  S)  shown  as  a  function  of  salinity 
S  with  the  Brillouin  frequency  shift  ub  as  a 
parameter  between  7.4  and  7.5  GHz  in  steps  of  50 
MHz. 
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Fig.  2.  Temperature  T(a)B,  S)  shown  as  a  function  of  the 
Brillouin  frequency  shift  ub  for  values  of  salinity  S 
from  30  96o  to  40  %o  in  steps  of  2.5  %o. 


Sound  velocity  [m/s] 
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Fig.  3.  Sound  velocity,  vs(vb,  S),  shown  as  a  function  of 
salinity  S  v>4th  the  Brillouin  frequency  shift  vb  as  a 
parameter  (in  steps  of  50  MHz). 
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Sound  velocity  [m/s] 


Fig.  4.  Sound  velocity,  vs(ub,  S),  shown  as  a  function  of  the 
Brillouin  frequency  shift  ub  for  values  of  salinity  S 
from  30  96o  to  40  %o  in  steps  of  5  %o, 
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Year  of  measurement 


Fig.  5.  VariabiUty  of  the  salinity  at  32n5'  S  17“30'  E,  in  the 
South  Atlantic,  near  the  Cape  of  Good  Hope,  at  a 
depth  of  10  meters.  Sampling  area:  30x100  km^. 
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Cumulative  fraction  of  locations 


Depth:  10  meters 


Fig.  6a.  Cumulative  histogram  of  the  variability  of  the 
salinity  standard  deviation  for  113  locations  for 
a  depth  of  10  meters. 
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Cumulative  fraction  of  locations 


Depth:  100  meters 
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Fig.  6b.  Cumulative  histogram  of  the  variability  of  the 
salinity  standard  deviation  for  113  locations  for  a 
depth  of  100  meters. 
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Fig.  7.  Seasonal  variability  of  the  salinity  at  34*00'  N  164*00'  E,  in 
the  middle  of  the  North  Pacific,  at  a  depths  of  10  m,  50  m, 
and  100  m.  Sampling  area:  50x50  km^.  The  standard 
deviation  represents  the  overall  variability,  the  seasonal 
trend  is  not  taken  into  account  in  the  calculation. 


On  Brillouin  UDAR ... . 


32/36 


Fry,  et.  al. 


Fig.  8.  The  uncertainty,  5T,  in  the  determination  of 
temperature  as  a  function  of  the  salinity  for 
5S»0.5%o  and  5vb=1  and  various  Brillouin 
shifts  from  7.4  to  7.7  GHz  in  steps  of  50  MHz. 
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Fig.  9.  On  the  left-side  y-axis,  the  uncertainty,  5T,  in  the 
determination  of  temperature  is  plotted  as  a 
function  of  the  salinity  standard  deviation,  for 
5ub=1  MHz  (solid  lines)  and  5ub=4  MHz  (dotted 
lines).  In  both  cases  we  have  taken  ub=7.5GHz 
and  S=35%o.  It  is  superposed  on  the  salinity 
standard  deviation  cumulative  probability 
distribution  at  100  meters,  on  the  right-side  y- 
axis  (same  as  Fig.  6b). 
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Cumulative  fraction  of  locations 


Sound  velocity  standard  deviation  [ro/s] 


Fig.  10  The  uncertainty  5vs  in  the  determination  of 
sound  speed  as  a  function  of  salinity  with  ub  as  a 
parameter  that  is  varied  from  7.4  GHz  to  7.7  GHz 
in  steps  of  50  MHz.  For  these  calculations  we 
have  talten  8ub*1  MHz  and  2B=0.5%o. 
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Fig.  11  On  the  left-side  y-axis,  the  uncertainty,  5vs,  in 
the  determination  of  sound  velocity  is  plotted  as 
a  function  of  the  salinity  standard  deviation,  for 
6'Ub=1  MHz  (solid  lines)  and  5ub=4  MHz  (dotted 
lines).  In  both  cases  we  have  taken  ub=7.5GHz 
and  S=35%o.  It  is  superposed  on  the  salinity 
standard  deviation  cumulative  probability 
distribution  at  100  meters,  on  the  right-side  y- 
axis  (same  as  Fig.  6b). 
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